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Abstract 

Low-capacitance Josephson junction arrays in the parameter range where 
single charges can be controlled are suggested as possible physical realiza- 
tions of the elements which have been considered in the context of quantum 
computers. We discuss single and multiple quantum bit systems. The sys- 
tems are controlled by applied gate voltages, which also allow the necessary 
manipulation of the quantum states. We estimate that the phase coherence 
time is sufficiently long for experimental demonstration of the principles of 
quantum computation. 

I. INTRODUCTION 

The issue of quantum computation has attracted much attention recently JlJ . Quantum 
algorithms can perform certain types of calculations much faster than classical computers . 
The basic concepts of quantum computation are quantum operations (gates) on quantum 
bits (qubits) and registers, i. e. arrays of qubits. A qubit can be any two-level system which 
can be prepared in arbitrary superpositions of its two eigenstates. Quantum computation 
requires "quantum state engineering" , i. e. the processes of preparation and manipulation of 
these quantum states in a controllable way. For quantum registers also "entangled states" 
(like the EPR state of two spins) have to be constructed. This necessitates a coupling 
between different qubits, and one has to be able to construct many-qubits states (many- 
particle states) in a controllable way. Several physical systems have been proposed as qubits 



lUU]. A serious limitation is the requirement that the phase coherence time is sufficiently 
long to allow the coherent quantum manipulations. The theoretically and experimentally 
most advanced system appears to be the chain with trapped ions || . 

In this paper we propose an alternative system, composed of low-capacitance Josephson 
junctions. The coherent tunneling of Cooper pairs mixes different charge states. By con- 
trolling the gate voltages we can control the strength of the mixing. The physics of coherent 
Cooper-pair tunneling in this system has been established before f5|-|7|. The algorithms of 
quantum computation introduce new, well-defined rules. Their realization in experiments 
creates a new challenge. Here, we consider a one-bit system and we describe the possible 
ways of constructing quantum states. Then we focus on a two-bit system, where we pro- 
pose a controllable coupling and we discuss the construction of two-bit states. In Section 
2, we consider ideal systems. In Section 3 we include the coupling to a realistic external 
electro dynamic environment which limits the phase coherence time. 

II. IDEALIZED ONE AND TWO-BIT SYSTEMS 

A. One-bit system and gates. 

The system we propose as a qubit is shown in Fig. |TJ. It consists of two small super- 
conducting grains connected by a tunnel junction. The junction is characterized by its 
capacitance Cj and the Josephson energy Ej. An ideal voltage source is connected to the 
system via two external capacitors C. We assume the total number of electrons in both 
grains to be even, so that at low enough temperatures all electrons are paired. Moreover, 
we assume that the superconducting energy gap A is the largest energy in the problem. 
Thus we can neglect pair-breaking processes. The complete set of possible quantum states 
is, then, characterized by and - the numbers of extra Cooper pairs on the left and 
the right island, respectively. Since there is no tunneling through the external capacitors, 
the total number of the extra Cooper pairs N = ul + ur is constant. In the following we 
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assume N = 0. Hence the set of basis states is parameterized by the number of Cooper pairs 
on one island or their difference n = (til — nn)/2. The Hamiltonian of this system is 

TT [2en-CV/2f 6 , 

H ~ C + 2C, - E "^*^>' (1) 
where G is the conjugate to the variable 2en || and $ — ~ is the flux quantum. We 
consider systems where the charging energy of the internal capacitor Ec 3 = (2e) 2 /2Cj is 
much larger than Ej. In this regime the energy levels are determined by the charging part 
of the Hamiltonian (|l]) for most values of the external voltage V. However, when V is such 
that the charging energies of two neighboring states |n > and \n + 1 > are close to each 
other, the Josephson tunneling becomes relevant. Then the two eigenstates are symmetric 
and antisymmetric combinations of |n > and \n + 1 > with an energy gap of Ej between 
them. 

We concentrate on the voltage interval near one such degeneracy point. We demand, 
however, the interval to be wide enough, so that for most values of the voltage away from 
degeneracy to a reasonable approximation the eigenstates are \n > and \n + 1 >. Near the 
degeneracy points we can translate (jl]) to the spin-| language: 

where | |>= \n > and | |>= |n + 1 >. Using this language we can immediately propose a 
few one-bit operations. If, for example, one chooses the operating point (i. e. the voltage) 
away from the degeneracy, the logical states |0 > and |1 > are just | |> and | f> respectively. 
Then, switching the system suddenly to the degeneracy point for a time At and switching 
suddenly back, we can perform one of the basic one-bit operations (gates) - a spin flip: 



t cos(^) zsin(^ 
V zsin(^) cos(^) j 



(3) 



To get rid of the time-dependent phases we perform all the calculations is the interaction 
picture, when the zero order Hamiltonian is the Hamiltonian at the operation point. To 
estimate the time-width At of the voltage pulse needed for a total spin flip, we note that a 



reasonable experimental value of Ej is of the order of IK. It cannot be chosen much smaller, 
since the condition k B T <C Ej must be satisfied. Therefore the time-width is very short: 
At w 10" 10 s. 

However there exists an alternative way to perform a coherent spin flip, as follows: The 
system is pushed adiabatically to the degeneracy point. Then the ac voltage with frequency 
Ej/h is applied. The process is completely analogous to the paramagnetic resonance (here 
the constant magnetic field component is in the x-direction, while the oscillating one is in 
the z-direction). The time- width of the ac pulse needed for the total spin flip depends on 
its amplitude, therefore it can be chosen much longer then lCT 10 s. 

One can also perform an operation which changes only the relative phase between the 
the states | |> and | |>. To this end one should suddenly switch to another voltage away 
from the degeneracy and switch suddenly to the original point. The phase is acquired due 
to difference in energy gaps between the two points. 

B. Two-bit system. 

To perform two-bit operations (gates) which result in entangled states, one has to couple 
the qubits in a controllable way. The best situation would be one where the coupling can 
be switched off, leaving the qubits uncoupled when no operations are performed. Unfortu- 
nately, this option appears difficult to realize in microscopic and mesoscopic systems. As 
an alternative we suggest for a system with a weak constant coupling between the qubits. 
We will show that by tuning the energy gaps of the individual qubits we can change the 
effective strength of the coupling. We propose to couple two qubits using an inductance as 
shown in Fig.(0). One can easily see that if L = the system reduces to two uncoupled 
qubits, while if L = oo the Coulomb interaction couples both qubits strongly. Therefore, 
we are interested in low values of L (the exact criterion to be specified). The Hamiltonian 
describing this system is 
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tt _ f (2en, - ViC t ) 2 . 9* J g 2 2 

e(ni - n 2 )q e 2 C t 



Cj ^K + ^r, (4) 

where g denotes the total charge on the external capacitors of both qubits, is its conjugate 
variable, and C^ 1 = Cj l + 2C~ 1 . The oscillator described by q and variables produces 
an effective mean field coupling between the qubits for frequencies smaller than ujlc = 



1 / y/2C t L. In order to have this coupling in a wide enough voltage interval around the 
degeneracy point, we demand that 

•4^»1. (5) 

Although the integration over q and can be performed exactly starting with (M), the 
resulting effective action is not simple. We prefer an alternative way. First we perform a 
canonical transformation 



q = q- 2e ^-"^ , = 0; 

hi = rii, Q l = Q l - {-!)* g< 



(6) 



i — 1, 2, which results in a new Hamiltonian (we omit the tildes): 



[ (2en t - CV.J2Y 



277 r e l + (-iy Ct 



%/2 V ' C 



.1 



2(2C t ) 2L 
Now we make an additional assumption 

<^>«(|t) 2 - w 

This assumption is, actually, necessary if one wishes to have Josephson tunneling terms in 
the Hamiltonian. Qualitatively, one can observe from (0), that if (|8|) does not hold, the 
Ej cos(..) terms are washed out. Below we will obtain this result in a more rigorous way. 

Assuming the condition (^) we expand the £7jcos(..) terms of (0) in powers of and 
we neglect powers higher than linear. Then we can perform the integrations over q and 0, 



which are again Gaussian. As a result we obtain an effective Hamiltonian, which consists of 
two one-bit Hamiltonians (|1|) and a coupling term: 

/ 2tt 



hit 



Ei 



sm 



$ /2 



61 — sin 



2tt 
%J2 



0, 



where 



E, 



•57T 



2 C t 2 EjL 



(9) 



(10) 



This Hamiltonian is valid for frequencies which are smaller than ujlc- I n the spin-| language 
it is 



tf„„ = -fK I, -0 2 



(ii) 



The mixed term in (|TT|) is important in certain situations. Assume that El is comparable 
to E] and the voltages V\ and V2 are such that both qubits are out of degeneracy. Then, 
to a good approximation, the eigenstates of the two-bit system without coupling are | ||>, 
I IT>> I Tl> an d I TT>- In a general situation, these states are separated by energies which 
are larger or much larger than Ej. Therefore, the effect of the coupling is small. If, however, 
a pair of these state is degenerate, the coupling may lift the degeneracy, changing the 
eigenstates drastically. For example, if V\ = V2, the states | ||> and | f|> are degenerate. 
In this case the correct eigenstates are: ||> +| tl>) an d 7^(1 !t> — I Tl>) with the 
energy splitting El between them. 

Now we are ready to propose a way to perform two-bit operations which will result in 
entangled states. For this we choose the operating points for the qubits at different voltages. 
Then we switch suddenly the voltages to be equal for a time At and switch suddenly back. 
The result is a "generalized" spin-flip, which may be described in the basis {| ||>, | IT>> 
I t!>i I TT>} by a matrix: 















cos(^) i sinf^) 



V 2h J 



Ozsin(^) cos(^) 
1 



(12) 
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The gate flTj ) together with the one-bit gates (|3]) constitute a universal set, i.e., they are 
sufficient for quantum computation. 

Instead of applying very short voltage pulses, one again can move the system adiabatically 
to the degeneracy point (Vi = V2), and, then, apply an ac voltage pulse in the antisymmetric 
channel V\ — V 2 = Aexp(iE L t/h). 



III. CIRCUIT EFFECTS AND DISSIPATION 

A. One-bit system. Energy relaxation out of the degeneracy point. 

The idealized picture outlined above has to be corrected with respect to the possible 
dissipation mechanisms, which, as usual, cause decoherence and energy relaxation. In this 
paper we focus on the effect of ohmic dissipation in the circuit, which originates mostly from 
the voltage sources. We also consider the effect of LC resonances in the circuit. We neglect 
quasi-particle tunneling since it is strongly suppressed at low temperatures and low voltages. 
The system is shown in Fig. |3|. We include the inductance L explicitly for two reasons. First, 
we would like to couple the bath of oscillators to a continuous charge variable, rather than 
directly coupling them to the discrete variable n. Second, we would like to treat the LC 
oscillatory mode separately from the bath oscillators, since it plays an important role in the 
two-bit system. The Hamiltonian for this system is: 

TT (2en - VC t ) 2 „ ( 6 \ q 2 2 2enq TT . . 



where 



#bath(g) = 



A. + ?M L - -A_ ? V 

2m j 2 \ 3 mjuy) 



(14) 



A 2 

and J(w) = I J2j ^^"^(^ — w i) = ^ w ■ m (HI) an d (0) ^ ne variable q denotes the charge 



on the external capacitors, while is its conjugate variable. The coupling to the bath fll4"D 
produces the usual dissipative term - Rq - in the equation of motion for q. 
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Although the Hamiltonian fll3|) is quite transparent, another equivalent form is more 
convenient. It is obtained after two consequent canonical transformations: 



q = q 



+ J2 j 



(15) 



^3 III . :C 



-q 



Pj = pj 



and 



q = q- 



n = n 



2enC t 
Cj 



e = e + g ( 



(16) 



Now the Hamiltonian takes the form: 



(2en - CV/2f 
H = — — Ej cos 



2n (n Ct ^ 



where 



^bath(0) = E 



3 L 



r P) , m J U j ~2 
— - — -X; 

2m, 2 3 



(17) 



2L 



(18) 



The two forms (|13D and (|17[) differ mostly in the definition of the phase variable conjugated 
to q. If L is small enough, the variable <fi fluctuates in a small interval around zero. However, 
the fluctuations of <f> are around a bath-dependent value, which may fluctuate strongly itself. 
Below we will see that both forms may be useful in different situations. 

First, we estimate the energy relaxation time, r r , due to the ohmic dissipation. We 
assume that the system is away from the degeneracy point and is prepared in one of its 
eigenstates (\n > or |n + 1 >) at the beginning. We use the Hamiltonian in the second form 
(|17|). The part of ([[7]) connecting the states |n > and |n + 1 > is 

. 2vr C t ~; 



Ht = Y exp (^®) exp 



+ h.c. 



(19) 



$o/2 Cj 

At this point we can directly apply the standard Golden Rule results for the transition rate 
from the state |n > to the state \n + 1 > p|-[T2| : 



7T 



(20) 



P(AE) 



K(t) 



2nh 



dt exp 
dw ReZ t (w 



(~i2 



10 id 

Z7 1 = iuC t + 



R K 
1 



c f(t) 
coth 



h 

Two 



AEt 



2k B T 



[cos(ut) — 1] — % sin(u;£) 



(21) 
(22) 
(23) 



R + iojL 

Here AE is the energy gap between the two states. The qualitative behavior of the system is, 
as usual, controlled by the dimensionless conductance g = j^. In our system a remarkable 

renormalization of the controlling parameter g occurs: from (pl|) one can observe that g = 

c 2 

^ig is the relevant parameter. Thus, making the external capacitances C much smaller than 
the internal one - Cj, we can drastically reduce the effect of the dissipation. Physically, this 
means that the fluctuations produced by the resistor are screened by the small capacitors, 
and, therefore, they have little effect on the Josephson junction. 

To be more concrete, we exploit the asymptotic formula for P(AE) fT2f 



P(AE) 



exp(-27/<7) 1 



7T AE 



gE 



2/a 



T(2/g) AE 

where T(..) is the Gamma function. For large values of g we obtain: 

1 g AE 



(24) 



T r = 



r(AE) 



T, 



op 



2vr 2 E 3 



(25) 



where r n 



is the operation time (see @). 



B. One-bit system. Circuit effects at the degeneracy point. 

At the degeneracy point the system is equivalent to the two-level model with a weak 



ohmic dissipation, which was extensively studied during the last decades (13|. It is well 
known that when g ^> 1 coherent oscillations take place. These oscillations make the spin- 
flip operation (|3[) possible. The decay time of the coherent oscillations is given by 

~ JL — - J- 

Td ~ 2vr 2 Ej ~ 2n* T ° p ' 1 ^ 

and the energy gap Ej is slightly renormalized: 



E^EA^)-' . (27) 



The physical cut-off u c is usually a system-dependent property. For a pure ohmic dissipation 
caused by a metallic resistor it may be as high as the Drude frequency. However, when 
additional capacitances and inductances are present in the circuit, the relevant cut-off may 
be lowered to the characteristic LC frequencies. 

As was said above, the LC phase fluctuations can wash out the Josephson coupling. To 
see this, we begin with the first form of the Hamiltonian fll3l ) and we integrate out the bath 
variables and the oscillatory mode variables - 0, q. The partition function reads: 



Z = J2 DnDQ x 

J no 



exp 
where 



no 

h 



/ dr(2eiQn- H (n,e))- / drdr' -G(r - T')n(r)n{T') 
Jo Jo Jo 2 



(28) 



and 

r( v _ 4C t e 2 1 

G[Un) ~ (l + C t Lu* + C t R\u n \) ■ m 

Below we will show, that in the relevant parameters' range, the following inequality holds 

I I R 

> -7= > - . (31) 



C t R y/LC t L 

Therefore, the natural cutoff for ( j30[) would be ui c = ujlc = v ^r - To make a rough 
estimation we replace ( |30D by 

G(u n ) = - 4 ^ 2 6 (1 - C t Lu) 2 n - C t R) , uj n <uj c 

G(uj n ) = , u n > u c . (32) 

The kernel fl3"2] ) is a sum of a trivial charging energy renormalization (the first term), 
the usual ohmic dissipation (the third term) and the inductive term, which we will 
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focus on. We apply the standard charge representation technique 0], expanding 
exp drEj cos(2tt 1 ^2) in powers of Ej and integrating over term by term. One 

obtains, then, a path integral over integer charge paths with instantaneous "jumps" be- 
tween the different values of n. Each such "jump" contributes a multiplicative factor of 
Ej/2h to the weight of the path. It turns out that the inductive term contributes another 
multiplicative factor for each "jump" , so that Ej is renormalized as: 

* - e > ^ • (33) 

One can immediately observe that the condition that Ej is not renormalized to zero coin- 
cides with the small fluctuations condition @. We would like to emphasize that the phase 
fluctuations which may wash out the Josephson coupling in (^) are related to the "weakly 
fluctuating" phase 0, rather than to the "strongly fluctuating" (see fll5D). Thus the effects 
of the inductance and the dissipation may be well separated in this regime. Another way 
of viewing this separation is by noting that the LC phase fluctuations are fast, therefore 
they effectively wash out the slower processes (like Josephson tunneling). On the other hand 
the phase fluctuations, caused by the resistor, are large only at low frequencies, thus they 
cannot wash out the faster processes. 



C. Two-bit system. Circuit effects. 

In this subsection we show that the coupling term @ survives when the ohmic dissipation 
is present in the circuit. The analyzed system is shown in Fig.@. We include the auxiliary 
inductances L a in order to couple the dissipation bathes to continuous charge variable, rather 
than directly to the principal variables n and m. Finally we take L a — > 0. The Hamiltonian 
is given by: 

^ f (2en t - V t C t ) 2 ( f \ qf 2en t q \ 

ff = L ^ E J cos 27r ^T7^ + 7^r ~ —F^- + #bath(%) 



i=l,2 



2Cj V $o/2; 2C t Cj 



, (01 + 02) 2 (01 - 02) 2 f , 

4L a + 4(2L + L a )' 1 ] 
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Here q\ is the charge on the left-most external capacitor, while g 2 is the charge on the right- 
most external capacitor. <f>i and 2 are their conjugate phases (weakly fluctuating ones). 
We make a canonical transformation 



ft = ft - ^ , & = & ; 

fii = ni , 0i = 0j + , 

(i = 1, 2) and the new Hamiltonian looks like (we omit the tildes): 



(35) 



[ (2en t - CV^/2) 2 

F = £ 2 H^ EjCOS 



2 - ^a-a 



$o/2 V C 



.1 



+ + ^bath (ft + ^T 2en i 



+ 2 ) 2 , (0!-(/. 2 ) 2 



(36) 



4L a 4(2L + L a ) ' 

Although the canonical transformation (35) looks equivalent to ([IB]), there is one essential 
difference. In fll6|) we used the phase variable which included the bath fluctuations, while 
in ( pop we keep the weakly fluctuating phases <pi and 2 . Next, we expand the _Ejcos(..) 
terms in <pi and 2 and we integrate over the baths and over the two oscillatory modes. 
Finally we assume L a — ► 0. One can directly check that in the nondiagonal channel the 
mixed term of (|9[) emerges, accompanied with a bunch of other terms which are small at 
low frequency. Thus, the two-bit coupling, crucial for the quantum computation, is stable 
under the influence of the dissipation. 



IV. DISCUSSION 

Several conditions were assumed in this Letter in order to guarantee the validity of 
our results. Here we state these conditions once more and discuss the appropriate ranges 
of parameters. We start with Ej ~ 1 K as a suitable experimental condition. To satisfy 
Ej Eqj we take Cj ~ 1CT 16 F, which is an experimentally accessible value. As we would 
like A to be large (|5|), it seems that L and C t should be as small as possible. However, the 
two-bit coupling energy, El fllPP, should be larger than the temperature of the experiment. 
Assuming a reasonable working temperature of 20 mK, we demand El ~ 0.1K. From ([S]) 
and ( |i0|) we get C t = ElC]A 2 /e 2 . To have a wide enough operation voltage interval we take 
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A « 10, and obtain C t ~ 10~ 17 — 10~ 16 F and L ps 10~ 8 — 10~ 7 H. Thus the renormalization of 
g is of the order of 10, and T r /r op w 10 2 — 10 3 (assuming the realistic value R ~ 100 O) (p5|). 
Finally we observe that in this range of parameters the inequalities (pj) and fl3~T|) are always 
satisfied. We conclude that the quantum manipulations we have discussed in this paper can 
be tested experimentally using the currently available lithographic and cryogenic techniques. 
Application of the Josephson junction system as an element of a quantum computer is a 
more subtle issue, demanding either the fabrication of junctions with Cj < 10~ 16 F, or a 
further reduction of the working temperature. 
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FIG. 3. A qubit with dissipation. 
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FIG. 4. Two-bit system with dissipation. 
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